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It is shown that simultaneously changing the scattering length of an elongated, harmonically 
trapped Bose-Einstein condensate from positive to negative and inverting the axial portion of the 
trap, so that it becomes expulsive, results in a train of self-coherent solitonic pulses. Each pulse 
is itself a non-dispersive attractive Bose-Einstein condensate that rapidly self-cools. The axial trap 
functions as a waveguide. The solitons can be made robustly stable with the right choice of trap 
geometry, number of atoms, and interaction strength. Theoretical and numerical evidence suggests 
that such a pulsed atomic soliton laser can be made in present experiments. 
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I. INTRODUCTION 

Solitons have applications in a wide variety of physi- 
cal contexts, ranging from water waves to photonic crys- 
tals P, 0, 01 For example, they have been used in 
transatlantic communications systems, where the need 
for expensive amplifiers mid- line in fiber optic cables that 
run over long distances is reduced or eliminated. By 
virtue of their many uses, as well as their mathemati- 
cal beauty, solitons are a continuing subject of vigorous 
research (see, for example, 0, ||)- In particular, they 
have proven highly useful in laser applications 0] . Non- 
linear materials are used to cause high-intensity coherent 
light waves emitted by lasers to self-focus into stable non- 
dispersive pulses. 

Bosc-Einstcin condensates (BEC's) arc coherent mat- 
ter waves in analogy to coherent light waves. BEC's are 
usually generated as a standing wave in a trap which 
functions as a cavity. When outcoupled from the trap, 
a BEC can provide a highly brilliant source of coherent 
matter-wave radiation, and as such is commonly called 
an atom laser. The challenge in making a BEC into a 
useful atom laser is in the outcoupling ((J Q. To this 
end, many experimental methods have been developed. 
Anderson and Kasevich @ tilted a BEC trapped in a pe- 
riodic potential created by a standing light wave. The 
gravitational field induced by the tilt caused the conden- 
sate to tunnel through the wells and interfere coherently, 
thereby creating a pulsed atom laser. Bloch ct al. used 
an external laser to change the spin state of atoms in two 
locations in a harmonic magnetic trap. The condensed 
atoms then spilled out, again due to gravitational effects; 
the two separate outcouplings allowed them to make the 
first clear demonstration of coherence along the whole 
length of the beam. Many experiments and proposals 
have followed. All of these atom lasers were studied in 
the context of repulsive BEC's. Repulsive BEC's natu- 
rally disperse in all directions; even axial confinement in 



a waveguide cannot prevent spreading in the direction of 
propagation. 

The unique contribution of the present study is to show 
how an attractive BEC can be made into an atom laser. 
A repulsive BEC fractures near conducting surfaces > 
spreads out, and in general is easily excited. In contrast, 
an attractive BEC, so long as it is axially confined and 
the experimental parameters are chosen properly, can be 
made into a pulsed atomic soliton laser which is robustly 
stable against all of these effects. Moreover, attractive 
BEC's in this form may be superior to repulsive BEC's 
in applications to atom chips [Til IT^ | and interferome- 
try [l3| . For instance, in the "nevatron", a BEC storage 
ring, the wave packet of repulsive Bose-condensed atoms 
circulates a few times before spreading out and decoher- 
ing This effect is accentuated by superconducting 
wires which lie transverse to the direction of propagation. 
A bright soliton would not only be non-dispersive, but, 
even if excited by the passage over the jump in potential 
created by the wire, would quickly self-cool by emitting a 
small fraction of its atoms, typically less than a fraction 
of a percent [Tol ITf3 | . It could therefore circulate indef- 
initely, subject to three-body recombination rates and 
other effects beyond those of the mean field 0] ■ 

In the following, we explain how to create such a pulsed 
atomic soliton laser. Bright matter- wave solitons have 
been created, both singly and in trains 0. It is 
shown that a combination of the experimental techniques 
of Refs. 0] and together with the right choice of 
parameters, suffices to create a pulsed atomic soliton 
laser from an attractive BEC. After presenting the basic 
method in Sec. [HI we illustrate its viability via three- 
dimensional simulations in Sec. IIIII Then, in Sees. IIVI 
and [yj the stability criteria and important dynamical 
features of the pulsed atomic soliton laser are explained 
in detail. In Sec. IVII the simulations are discussed in 
light of the results of Sees. IIVI andlVl Finally, in Sec.lVlTl 
we conclude. 
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II. BASIC METHOD 

The 3D Nonlinear Schrodingcr equation (NLS) or 
Gross-Pitaevskii equation which describes the mean field 
of the BEC is written as [l7| 



2m 



V 2 +gN\V\ 2 + V(r) 



where 



(1) 



(2) 



g = 4nh 2 a/m, a is the s-wave scattering length, m is 
the atomic mass, N is the number of condensed atoms, 
the condensate order parameter <!> — ^(r,t) has been 
normalized to one, and axisymmetric harmonic confine- 
ment has been assumed. Note that for negative scat- 
tering length, or attractive nonlinearity, solutions are li- 
able to collapse in certain parameter regimes 0, EH , as 
shall be discussed below. With the exception of Sec. IV Bl 
where the decoherence time between pulses is estimated, 
it is assumed that the BEC is described by Eq. Q. 

The basic method for creating the pulsed atomic soli- 
ton laser is as follows. 

1. A repulsive BEC is created in an elongated har- 
monic trap such that uj z <C u> p . 

2. The scattering length is made small and nega- 
tive via the now well-established experimental tech- 
nique of using a magnetically induced Feshbach res- 
onance m Him nil. 

3. Simultaneously, the axial potential is changed from 
small and attractive (u> 7 , re al) to small and expul- 
sive (u> z imaginary) [n^ lisl I23] . 

4. The condensate becomes modulationally unstable 
to spatial pulse formation. This instability is 
non-dissipative. The pulses are seeded by self- 
interference of the order parameter, as we have else- 
where described [24| • The initial growth rate of the 
modulational instability can be calculated via lin- 
ear perturbation theory |25j . 

5. The ensuing solitonic pulses are subject to pri- 
mary collapse in two or three dimensions, as well 
as secondary collapse due to soliton-soliton inter- 
actions |24|. Furthermore, if they are too long 
in the ^-direction they can be torn apart by the 
force of the expulsive harmonic potential [l^ . With 
the right choice of parameters, these effects can be 
avoided and the solitonic pulses made robustly sta- 
ble, as shall be described in Sec. IIVI 

6. The solitons continue to accelerate. Their relative 
spacing increases as Az oc exp(|w 2 | t). 
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FIG. 1: Shown is the evolution of an attractive Bose-Einstein 
condensate into a pulsed atomic soliton laser. Time slices 
of the line density in z are shown for x, y = 0. Modula- 
tional instability of the initial density profile is seeded by self- 
interference of the order parameter, so that solitons form first 
at the cloud edges and later towards the center. The lat- 
est, top panel, shows that a well-separated set of stable soli- 
tonic pulses are produced. Note that, for N = 10 4 atoms, 
a — —3ao, and a trap geometry of uj p — 2ir x 2.44 kHz, 
uj z = 2ni x 2.26 Hz, the time units are scaled to 22 ms and 
the spatial units to 10 /im. 



7. The expulsive harmonic potential can eventually be 
coupled to a linear or even a flat potential for appli- 
cations. We do not discuss the various possibilities 
here. 

An important point is that even if the coupling is not 
smooth, each soliton responds to perturbation by a shift 
in its phase and by emitting a small number of atoms, 
typically a fraction of a percent of the total number in 
the soliton. Insofar as an excited BEC described by the 
NLS models a condensate plus thermal field [2ti |27| . 
where the "condensate" is a stationary solution, one may 
term this process self-cooling. As each soliton is itself 
a BEC, this model can be applied to each pulse sepa- 
rately. Self-cooling to T = in an expulsive harmonic 
potential occurs exponentially, with the density fluctua- 
tions at the center of each solitonic pulse falling off as 
exp(-K|t) 0. 
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FIG. 2: Shown are the evolution of the density and phase along a two dimensional cut at y = for the simulation of Fig. Q 
A set of well-defined solitonic pulses is evident in the latest (top) panel. The strong variations in the phase at late times is 
due to the high momentum of the solitons caused by the expulsive harmonic potential. Note that the phase is shown on the 
color circle, i.e., modulo 2ir, while the density is in arbitrary relative units rescaled for each plot. The aspect ratio of the plots 
showing a region of 0.822 by 153 length units was changed for visualization; length and time units are the same as in Fig. 



One may ask what this schemata has in common with 
the operation of a normal, light-wave laser. Insofar as the 
initial condition is an excited mode of the harmonic trap 
plus mean field potential, and the ensuing pulse train is 
a much lower energy mode, this situation has a certain 
analogy with population inversion. As was mentioned 
in the introduction, the initial trapping potential may 
be considered as a cavity, with the outcoupling provided 
by the sudden change to an expulsive trapping poten- 
tial in the axial direction. However, the emission of soli- 
tonic pulses is not stimulated, as strictly required for the 
use of the acronym LASER (light amplification by stim- 
ulated emission of radiation), but rather spontaneous. 
The analogy of the proposed scheme to the operation of 
a laser is therefore rather in the output than in the de- 
tailed mechanism of its operation: one produces a train 
of non-phase-locked self-coherent pulses. 



III. PULSED ATOMIC SOLITON LASER 
DYNAMICS: PROOF OF PRINCIPLE 
SIMULATIONS 



Three dimensional simulations of Eq. Q were per- 
formed, with parameters which satisfied the criteria given 
in Sec. IIVI Cylindrical symmetry was assumed in or- 
der to make computations with a large grid size possi- 
ble (2048 x 16) [23. The initial profile was obtained by 
imaginary time relaxation. This resulted in a Thomas- 
Fermi-like profile in the z direction (see Eq. 112611 below 
and [13) > an d a nearly Gaussian one in x and y. A trap 
of aspect ratio lo p /uj z = (W^) 2 = 538 and nonlin- 
carity parameter a s oN/£ z o = 1.02 was used to produce 
the initial state, where w z o, £ z o an d a s o a ll refer to these 
initial conditions, and £i = y^U/ (mu>i). 

At t = 0, the longitudinal trapping frequency lo z q 
was changed to a weak expulsive harmonic potential 
with lu z = 0.5iuj z o, and the nonlinearity was switched 
from repulsive to attractive, with —aN/£ Z Q = 0.0854 = 
0.00368 £ z o/(p. As attractive BEC's can collapse in three 
dimensions, this is an important point in the choice 
of parameters. The length unit in the simulations is 
u = 0.56€ z q an d the time unit r = 2mu 2 /h. To com- 
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pare with experimental parameters, one must choose the 
number of particles and a scaling factor, e. g. u = lOfim 
and N = 10 4 , which corresponds to r = 22ms, uj z q = 
2?r x 4.53Hz, u z = 2iri x 2.26Hz, and uj p = 2tt x 2.44kHz 
with a scattering length of a = — 3ao- 

Figures Q an( l HI illustrate the evolution of the density 
and phase of the condensate in time slices through the x-z 
plane for y — with the above described initial state and 
parameters. Several observations may be made based 
on the figures. Firstly, the final number of solitons is 
14. Secondly, they are stable against collapse and, once 
formed, do not subsequently interact over the lifetime of 
the simulation. Thirdly, solitons form first at the edges 
of the cloud, then later towards the center, as was also 
observed in a simplified model in our previous work |24j | . 
We note that no white or colored noise was added to this 
simulation. The reasoning behind our choice of parame- 
ters will be discussed in Sec. IIVI while the details of the 
simulation itself will be interpreted in Sees. IVl and IVII 



IV. STABILITY CRITERIA 

In order that the pulsed atomic soliton laser be ro- 
bustly stable over the lifetime of an experiment, a set of 
criteria must be satisfied. These criteria are detailed in 
the following subsections. Note that in the below con- 
siderations wc arc interested in stability for experimental 
purposes, not mathematical stability to infinite time. 



A. Two-dimensional primary collapse 

Two-dimensional transverse primary collapse must be 
avoided. In the case of strongly anisotropic axisymmct- 
ric confinement, one may adiabatically separate the slow 
longitudinal from the fast transverse degrees of freedom. 
The adiabatically varying transverse state obeys a 2D 
NLS which shows an instability towards collapse. The 
criterion for stability found by numerical integration of 
the 2D NLS is 



8naniY)(z, t) < rj, 



2D 



11.7. 



(3) 



where n -mt z, t) is the local axial line density of the con- 
densate I23. If adiabaticity is violated, collapse can also 
happen at weaker nonlincarity due to transverse oscilla- 
tions on a time scale tt/uj p [3(j. When the longitudinal 
dynamics is significantly slower than this time scale, the 
adiabatic separation of scales is valid. If, additionally, the 
transverse nonlinearity is weak, i.e., \8Traniu(z,t)\ -C 1, 
the longitudinal equation reduces to the quasi-lD NLS 



■— d 2 z +g 1D N\i;\ 
2m 



1 22 
— mui,z 
2 



^ = ihdti>, (4) 



where giD = 2 au> n Ti is the renormalizcd 
pling constant [l(j, provided l p ^> \a\ > 
y/h/(muj p ). 



uasi-lD cou- 
with 



' p — 



B. Three-dimensional primary collapse 

Three-dimensional primary collapse of the individual 
solitonic pulses in the atom laser must be avoided. The 
static condition based on imaginary time relaxation of 
Eq. is 



< r)f = 0.627 . 



(5) 



where iV a i s ~ N/N s is the number of atoms in a soliton, 
with N s the number of solitons (see Sec. IVII below) . Re- 
placing the < sign with a -C sign ensures stability for 
an excited soliton. Note that the value of t^P can vary 
slightly as the anisotropy of the trap changes [liit l32j |. 
However, 3D collapse is essentially an isotropic phe- 
nomenon, with the trap simply setting the initial con- 
ditions l33l. 



C. Explosion of individual solitonic pulses 

The soliton can become unstable when the expulsive 
potential overcomes the balance between the mean field 
energy and kinetic energy necessary for the soliton's ex- 
istence and destroys it by tearing it in two [l(| . We term 
this kind of possible instability explosion. In order to 
avoid explosion, the geometry must be chosen so that 
£ z ^S> 4oi, where l z = y/n/m\uj z \. The soliton length 
4oi ~ where £ is the healing length, can be deter- 
mined as follows. Taking the form of the soliton as the 
well-known solution in one dimension for a constant po- 
tential (see and references therein) 



1 



: sech 



(6) 



and substituting into Eq. (0J while temporarily ne- 
glecting the trapping potential, one obtains t z = 
2ft 2 /(m|£i D |iV a i s ). Then 



tsol 



\a\N K 



' p ■ 



(7) 



where the wavefunction has been renormalizcd to account 
for the division of TV total atoms into iV a i s atoms in any 
given soliton. Note that the prefactor is the inverse of 
the ratio which must be small to avoid 3D collapse, as 
given by Eq. JSJ in the preceeding subsection. 

A variational analysis based on a hyperbolic secant 
ansatz and Eq. Q with the trapping potential gives a 
more precise condition to avoid explosion as |l6j | 



> 



I .SOl 



2 6^4xl/4 



3 3 



= 2.20. 



(8) 



In the quasi- ID regime far from collapse, such a varia- 
tional analysis typically gives estimates to better than 
1%. As in Eq. JSJ, the > sign can be replaced with a ^ 
sign to ensure stability for an excited soliton. 
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D. Soliton— soliton interaction and secondary 
collapse 

The harmonic potential must be sufficiently strong so 
as to prevent secondary collapse caused by soliton- soliton 
interaction. If two solitons overlap coherently they can 
violate the stability criterion of Eq. due to the dou- 
bling of the number of atoms. Even if, due to decoherence 
during soliton propagation (see Sec. lVBl below). their rel- 
ative phase is not defined prior to interaction, upon in- 
teracting they develop a well-defined relative phase (34|. 

In order to treat soliton motion in a slowly varying po- 
tential, as defined explicitly by the condition of Eq. JHJ), 
one may suppose a separation of scales, as may be for- 
mally defined by a multiscalc analysis (see, for example, 
Ref. |35l|). In this case an approximate equation of mo- 
tion for the relative soliton parameters is given by |3^ . l3^ | 

1 8?iV • ( ^ 

"> = 77— sm(0) cxp(— nr ) 



ni- 



87Tk 



2 ,.3 



where 



cos(0) exp(— nr) — ui^r , (9) 



-Ng m m/h 2 , (10) 



4> = (f>i — 02 is the relative phase, r = \zi — Z2I is the rel- 
ative position, with the two solitons indicated by indices 
1,2, and motion according to Eq. has been assumed. 
Equation @ describes a separated soliton pair, i.e., the 
motion outside the interaction region: it breaks down as 
they overlap fully. With respect to Refs. j3(||33], we have 
here added the physical units relevant for the BEC and 
the expulsive harmonic potential. 

To prevent soliton-soliton interaction it is necessary 
that the potential due to the expulsive harmonic poten- 
tial be much stronger than that due to the attraction be- 
tween solitons. Taking cf>(t) = <fi(t) = 0, which assumes 
that the solitons are initially in-phasc and have the same 
amplitude, the two potentials are given by 



1 



2 2 



sol 



u 2 k 2 , x 

> 77- exp(-Kr) 

m 



(11) 
(12) 



In case the solitons are not initially in phase or do not 
have the same amplitude, the criterion will only be less 
stringent. It is therefore sufficient that 



t: 



aN a . 



AT? 



26 , 
> —77 cxp(- 



■8ttN s /N^ s ). (13) 



Here the first factor in parentheses is the trap aspect 
ratio while the second factor is again the inverse of the 
3D collapse criterion of Eq. (J3J) ■ 

The question then arises as to whether or not the trap 
can be made sufficiently strong so as to prevent soliton 



interactions, as required by Eq. I|13|) . an d at the same 
time sufficiently weak so as not to cause the individual 
solitons to explode, as required by Eq. Putting these 
two criteria together, one finds 



X 2 > -« exp(-87rx) ; 



(14) 



where \ = N s /N a i s . This relation is always fulfilled, 
showing that the two criteria are compatible. 

V. DYNAMICAL FEATURES 

There are two aspects of the dynamics which are nec- 
essary to discuss in detail. Firstly, there has been some 
debate as to the mechanism of soliton formation. In 
Sec. IV Al it is argued that both dynamically generated 
fluctuations from self-interference of the order param- 
eter [24| and noise due to thermal fluctuations or 
fluctuations in the trappin g p otential , as seen exper- 
imentally close to surfaces [4t| , cause the BEC to become 
modulationally unstable on approximately the same time 
scale. Secondly, outside of the mean field model encap- 
sulated in the NLS, one may ask how long it takes for 
the relative phase of solitonic pulses in the atom laser to 
randomize, or decohere. In Sec. IV Bl an estimate of this 
time scale is made. 

Two issues which we do not discuss in any detail are 
quantum evaporation and center of mass motion. The 
former is studied in detail in Ref. |l6l |. where it is shown 
that matter- wave bright solitons in an expulsive potential 
evaporate and eventually explode. However, the tunnel- 
ing rates are so small in the parameter regimes of interest 
to the present work so as to be unimportant. With re- 
gards to the latter, in a harmonic potential the center 
of mass and relative degrees of freedom are entirely de- 
coupled, so that we need only consider the relative soli- 
ton motion |4l| . The center of mass motion is, in any 
case, trivial: z com (t) = z com (0) expd^lt) in the quasi- 
1D regime. 



Seeding of modulational instability: 
self-interference vs. noise 



In order to understand the mechanism of modulational 
instability for a non-uniform initial density profile and in 
the presence of a non-constant potential, it is necessary 
to briefly review modulational instability in the uniform 
case, which is well known from fiber optics [2f|. A linear 
response analysis reveals that, for attractive nonlincarity, 
a small sinusoidal modulation of a uniform state ipo with 
wavenumber k grows with time at a rate 7 given by 



T 



4m 2 



2to|c/id| "id 



™ 2 dI<?id| 2 



The maximum growth rate 



7, ng = 2u p \a\ nm 



(15) 



(16) 
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FIG. 3: Shown are the self-interference fringes of the or- 
der parameter which seed modulational instability, according 
to the Feynman propagator for a harmonic oscillator. The 
initial density profile was Thomas-Fermi; shown is the ratio 
of the density 66 ms later to the original density. Since the 
wavelength of the instability must be on the order of 2n£, 
where £ is the healing length, solitons form first on the edges 
of the cloud, due to the early long wavelength fringes in this 
region. At later times the wavelength of the fringes in the 
center also becomes longer. Shown is the linear equivalent 
of the panel depicting t — 3 in the full simulation of the 3D 
Gross-Pitaevskii equation illustrated in Figs. and [5] all pa- 
rameters are the same as the simulation, with length units 
scaled to 10 /im. 



the diffraction pattern of a rectangular initial density pro- 
file |43| . Our previous analysis was performed for a rect- 
angular initial density profile in order to obtain closed 
form analytic results |24j . The Feynman propagator is 
defined by 



tp(z,t)= dz'G(z,t;z',0)ip(z,0). 



(18) 



For a harmonic oscillator, the propagator is 



exp {i(z 2 - 2zz'/cost + z' 2 )/(2P z tanr)} 



i z \J 2iri | sin : 



where 



UJ z t . 



(20) 



In the limit t« 1, and for a rectangular initial density 
profile, the result of the integration of Eq. I|18|) can be 
Taylor expanded as 



hA(M)| 2 /l^,0)| 2 



1 



8l 2 T 



sin(fc + z + 5 — j ) sin(£;_z + 5 



L + 2z 
L 2 + 4z 2 



+ 



L-2z 
cos[(fc + — k-)z] 
(L + 2z) 2 (L-2z) 2 + (L + 2z)(L-2z) 



,(21) 



is obtained at wavenumber 



1/e, Z = £ p /y/4\a\n 1D . 



(17) 



where £ is the effective ID healing length of the con- 
densate an d n iD = ^IV'ol 2 = N/L is the line 
density. Growth occurs only if j 2 > 0, which implies 



< k < k n 



V2k„ 



This means that nonlinear fo- 



cusing can only be seeded by modulations of sufficiently 
long wavelength and is fastest at the length scale of 2tt£. 

For a non-uniform initial density profile, there are two 
ways in which modulational instability can occur. Either 
it can be seeded by noise, or it can be seeded by fringes 
caused by self-interference of the order parameter. The 
time scale of the two mechanisms turns out to be approx- 
imately the same, as shall be shown in the following. A 
similar number of solitons results, but in the former case 
they form first in the higher density regions (typically 
the center, for a Thomas- Fermi- like initial profile H3)i 
while in the latter case they form first at the edges, as 
illustrated in Figs. ^ and |21 and explained in our previous 
work 0. 

Consider first the case of self-interference. An anal- 
ysis based on the Feynman propagator for the linear 
Schrodinger equation in a harmonic potential shows that 
self-interference, or diffraction, of the order parameter 
leads to fringes which have the correct length scale to seed 
modulational instability j24[. Our findings of Ref. |24| 
were supported later by Kamchatnov et at, who used 
Whitham theory to describe the nonlinear evolution of 



sec(r)z ± L 
2l 2 sin(r) 



L 2 cot(r) 
87? 



\z\< 



L 



(22) 



To linear order in r, the trapping frequency drops out 
of the equations, since 1 2 t — Td/ra. Equations l|21() 
and (|22|) describe the formation of fringes. Note that, 
according to the argument of the exponential in the Feyn- 
man propagator l|19|) , at the quarter period the wavefunc- 
tion is fourier transformed with respect to its initial state. 
Therefore any initial wavefunction excepting a Gaussian 
must develop fringes. A time scale can be estimated from 
these prefactors in the expansion of Eq. (|21f) . Fringes ap- 
pear at a length scale £f r i n go at time 



t 



*- ''-fringe ■ 



(23) 



This argument can also be made simply by the units in 
the problem. The length scale at which modulational 
instability is maximally probable is ^fringe = 2irt;. There- 
fore, the time scale for fringe formation leading to mod- 
ulational instability may be estimated as 



rmgc 



2uj p \a\n 1I ) 



(24) 



where Uyd is the mean linear density and we have taken 
the mean density as n = fim/irl 2 in order to calcu- 
late the healing length. For the parameters of Sec. IIIII 
^interference — 41 ms. This is approximately the correct 
time scale, as observed in Figs. ^ and [21 
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In order to study the problem with a more realistic 
model than an initially rectangular density profile, the 
longitudinal variation of the density profile is taken as an 
inverted parabola. This is characteristic of the Thomas- 
Fermi limit in a harmonic trap, and is the generic experi- 
mental case 0| • At the same time, the transverse wave- 
function is taken as a Gaussian, in keeping with the quasi- 
1D approximation. The density then takes the form 



|*(f,0)| 2 
|^(z,0)| 2 
R 



IV^,o)| 5 



i 



■ exp 



x 2 + y 2 

' e 2 



t 2 {R 2 -z 2 ) 



U\\a\ 
(3N\a\£l 

\Jfijra\u); 



(25) 

(26) 

(27) 
(28) 



where N\ip(z)\ 2 is the longitudinal line density, R is the 
Thomas- Fermi radius, and £ z is the longitudinal oscilla- 
tor length. The linear development of the wavefunction 
may be found at any time by numerical integration of 
Eq. |JTSJ. Note that, in this case, u> z is imaginary for 
the expulsive harmonic potential. An example relevant 
to Sec. IIIII is shown in Fig. [3J The longer wavelength 
fringes are clearly visible near the edges of the cloud, as 
discussed in our previous work |24j . This leads to soli- 
ton formation near the edges of the cloud at early times 
and in the center at late times. The figure uses the same 
parameters as the simulations of Sec. IIIII and may be 
compared to the fourth panel from the bottom, or t = 3, 
in Figs. [I] and 

Consider now the case of modulational instability 
seeded by noise, rather than interference fringes. There 
are two kinds of noise. They originate in different phys- 
ical mechanisms. The first is classical white or coloured 
noise, which may be induced, for example, by fluctuations 
in the trapping potential. The second is thermal quan- 
tum noise, which corresponds to a thermal distribution of 
Boguliubov excitations. One may estimate the relevance 
of the latter from first principles. The Boguliubov quasi- 
particle dispersion relation for a Thomas-Fermi profile 
is 17] 



h 2 k 2 


\h 2 k 2 


J 2m 


_ 2m + 


h 2 k 2 


2 2 

muj v 


2m 


2 



2gN\^(r)f 



\r\ < R 
r\> R. 



(29) 

Substituting the wavenumber of maximum growth for 
modulational instability, Eq. I|17|) . into Eq. <|29l) . the re- 
sulting energy is 



mg 



V3.gr 



(30) 



where n is the mean density which can be estimated 
from Eq. I|25|) . An experimental situation may, e.g., 




aAaaIJJaaa a 






FIG. 4: Shown is the evolution of the density along a one- 
dimensional cut at x, y — 0, with the same parameters as 
the simulation of Fig. Q but with the addition of noise, as 
described in the text. The time scale is shorter than that 
observed in Fig. and the solitons form first in the center 
of the cloud, rather than on the outside, but the end result 
is the same: a set of well-defined solitonic pulses is evident 
in the latest (top) panel. The length and time units are the 
same as in Fig. 



correspond to an initial temperature of the condensate 
of T ~ Tbec/2, with Tbec = hij[N/C(3)} 1/3 and ui = 
(uipUizo) 1 / 3 . In this case, one can estimate the probabil- 
ity of a Boguliubov mode of the appropriate wavelength 
to seed modulational instability from the bosonic number 
distribution function 



i(E) = 



1 



cxp(E/k B T) + 1 



(31) 



Using the numbers from Sec. IIIII one obtains T ~ 0.14/iK 
and E^°s/k B ~ 0.23^K, so that n(E%*) ~ 0.17. Thus 
noise caused by Boguliubov fluctuations is present with 
a non-negligible occupation number for the parameters 
we have chosen. 

A time scale for the growth of seeding fluctuations can 
be estimated as 



t, 



1/7* 



1 



2uj p \a\n 



ID 



(32) 



This appears to be smaller by a factor of ir 2 than tf r i ng c • 
However, given the qualitative nature of the two esti- 
mates, and the fact that they have the same parameter 
dependence, it is not possible to state that noise defini- 
tively dominates over the fringe mechanism. To test this, 
we performed additional simulations with thermally dis- 
tributed and various other realizations of noise in the 
initial condition. These simulations show that the final 
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result is essentially the same as that shown in Figs. [3] 
and ^ 11 solitons result rather than 13, with forma- 
tion first in the regions of higher density and then lower, 
in contrast to the case of fringes alone. Figure (QJ il- 
lustrates an example of the early time evolution of the 
density. Solitons begin to form at t ~ 30 ms, in contrast 
to Fig. 2] where fringes appear to self-focus at t ~ 40 ms. 
We conclude that the two mechanisms do indeed coexist. 

Noise was added into the simulations in the follow- 
ing manner. Because the longitudinal and transverse de- 
grees of freedom are represented differently in the La- 
guerre DVR cylindrically symmetric algorithm we used, 
and most of the grid is not occupied in the initial state, 
one generates a great deal of high energy and high fre- 
quency oscillations with usual noise schemes, such as 
adding a small random number to the wavefunction or 
its Fourier transform. The fluctuations that would seed 
modulational instability are of long wavelength. There- 
fore, starting from the wavefunction on the grid in po- 
sition space an FFT was implemented in the longitudi- 
nal direction. Then the wavefunction was multiplied at 
each point by I + nr, where r was a random number be- 
tween —0.5 and 0.5 and n was the noise level. In Fig. 21 
r = 0.1 was used. In the transverse direction, the noise 
was added on only half the grid closest to the center. 
Finally, in order to allow the noise to "thermalize" as 
much as possible, the wavefunction was propagated in 
real time for positive scattering length, until the noise 
had fully distributed itself, i.e., for times much greater 
than 2it/\lu z \. This was intended to represent, qualita- 
tively, a semi-classical approximation to a thermal dis- 
tribution |2(| of Boguliubov modes. Figure 0] then 
follows the real time evolution starting with this initial 
wavefunction after the scattering length is turned nega- 
tive and the trap is changed to be longitudinally expul- 
sive. In simulations with smaller noise levels we observe 
a coexistence regime and crossover of both seeding mech- 
anisms, as the growth of thermal fluctuations is signifi- 
cantly delayed when they are initially very small. 



B. Phase decoherence time 



where t is the decoherence time. The chemical potential 
may be determined from Eq. i.e., in the quasi-lD 
approximation, to be 



/' 



-htOn 



N 

K 



1. 



(34) 



where it has been assumed that N s solitons of equal am- 
plitude are formed. Then, from the derivative of Eq. i|34|> 
with respect to N, 



Afi ^ 111 
~AN ~ iV 



(35) 



For Poissonian number fluctuations, one may take AN = 
\Z~N. Setting A(f> = 2tt, which is a measure of complete 
uncertainty in the relative phase and therefore decoher- 
ence, Eqs. l|3*3^) and (|3*5|l yield 



t = 



nhVN 



l/'l 



Substituting Eq. O into Eq. ®, 



21^ 

TV 



AM ^ — . 



(36) 



(37) 



We note that, in contrast to a repulsive condensate in 
Thomas- Fermi limit, for which the decoherence time 
t oc TV 1 / 10 , in the case of solitonic pulses formed by mod- 
ulational instability t cx TV -3 / 2 . However, unlike in the 
repulsive case, the number of atoms is limited by the 
collapse conditions of Sees. IIV Al and IIVBI For the 
parameters of Sec. IIIII the phase decoherence time may 
be calculated to be about 540 ms, so that the solitonic 
pulses shown in the figures are expected to be coherent 
over the evolution period depicted. 



VI. DISCUSSION OF SIMULATIONS: NUMBER 
OF SOLITONS AND REFINED STABILITY 
CONDITIONS 



A condensate adiabatically split into two halves on a 
time scale much shorter than the quantum revival time 
has an initially well-defined relative phase 0, . Esti- 
mates for the decoherence time 01 have been made in 
a number of specific contexts in the literature, as for ex- 
ample in the two- well problem |47l Esj or for two spin 
states [i^. A general discussion of this issue may be 
found in Ref. [34]]. Here, we follow the straightforward 
estimates made in a recent article on atom interferome- 
ters using Bose-Einstein condensates, in which the phase 
decoherence time was studied experimentally |5(J, lU| ■ 

The Schrodingcr phase of each soliton may be esti- 
mated from its wavefunction, which is proportional to 
exp(— i/Mt/h): 



A(j> = tAn/%, 



(33) 



In experiments, a good model of the initial state of the 
condensate when the scattering length is changed from 
positive to negative is a longitudinal Thomas-Fermi den- 
sity profile [13 ■ In the following, explicit estimates for 
the number of solitons generated by such a profile and 
criteria to avoid collapse, in terms of the parameters of 
a possible experiment, is compared to the more idealized 
situation discussed in Sec. IIVI 

Under the condition that a suitable seed for the modu- 
lational instability is provided, one can estimate the num- 
ber of solitons generated for an initially homogeneous 
profile along the z direction of length L by 



(38) 




7^ 
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where the modulational instability is assumed to take 
place at the wavelength of maximum growth given by 
Eq. 

The 2D collapse criterion for the initial state J3J can be 
refined by demanding that the solitons formed by modu- 
lational instability are themselves stable against 3D col- 
lapse and satisfy Eq. (jSJ). For simplicity, it is assumed 
that the initial condensate is split up into N s solitons 
of equal amplitude. As seen from the numerical simula- 
tions of Sccs. lllll and IV Al this is not strictly true, but it 
serves as a useful order of magnitude estimate. For the 
homogeneous initial profile, one finds from Eq. 138|) 



87r\a\N/L < -{^f 



1.0. 



(39) 



This estimate assumes a quasi-lD initial state, where 
the transverse trapping is tight, so that \a\ <C £ P <C £. 
Note that under these conditions the transverse oscillator 
length l p does not enter the collapse criteria for homoge- 
neous initial density profiles. 

In the case of an inhomogeneous initial density profile 
the above estimates can be generalized by assuming that 
the length scale of 27r£ for the modulational instability 
is still valid locally. The number of solitons can thus be 
estimated as 



N H = 



dz 



2tt£(z) 



(40) 



For the Thomas-Fermi density profile of Eq. one ob- 
tains 



criterion |yj for the initial wavefunction yields the same 
scaling but a prefactor which is an order of magnitude 
larger, i.e., a less stringent constraint. 

One may now compare these refined estimates to the 
simulations of Sees. IIIII and IIV Dl First, the choice of 
—aN/£ Z Q was taken to be about one order of magitude 
smaller than the upper limit given by the criteria to avoid 
collapse, according to Eq. l|l2")l. Secondly, Eq. (|4*T|l pre- 
dicts an upper bound on the number of solitons to be 
N s sa 13.3; in the simulations, between 11 and 14 solitons 
were observed, depending on the noise level. Thirdly, 
all forms of collapse have been successfully avoided, in- 
cluding soliton-soliton interactions which might lead to 
secondary collapse. 

Many other parameter regimes were studied numeri- 
cally. It was found that, for a rectangular initial profile 
and no noise, as was studied analytically in Ref. |24j . 
increasing the strength of the nonlinearity to the critical 
value of 77hom = 1 (sec Eq. I|43|) below) brought about im- 
mediate collapse at the borders of the condensate. That 
is, the first soliton formed collapsed. An order of mag- 
nitude decrease in ryhom to 0.1 led to delayed collapse 
which occurred after all solitons had been formed, while 
for 77hom = 0.01 no collapse occurred. Note that a rectan- 
gular initial density profile may be created by optically 
induced potentials which form end-caps , as were used 
in the experiment of Ref. |l9j | . 



VII. CONCLUSION 



N, 



TF 



3A|q| 4 

£ 2 



(41) 



Similarly, one can derive a collapse criterion for the 
Thomas-Fermi case. The aspect ratio of the trap enters 
explicitly, and the criterion to avoid collapse is given by 



lolJV 



(n 3D ) 3 £ £ 



These estimates restrict the number of particles for a 
given scattering length that can be used in a particular 
trap geometry. Combining the above results, one finds 
an upper bound for the number of solitons that can be 
generated from a given trap geometry used for preparing 
the initial state. In the case of a rectangular initial state, 
or homogeneous case, one obtains 



AT hom < 0.0635- 



L 



(43) 



For a condensate initially described by a Thomas-Fermi 
profile one finds 



Nj F < 0.055 ( ^ 



(44) 



The bounds l|43|) and (|44|) were based on the 3D collapse 
criterion of Eq. JSJ). A similar analysis based on the 2D 



We have shown both numerically and analytically that 
a pulsed atomic soliton laser is viable. In particular, 
the figures illustrate the evolution of such an atom laser 
with a set of realistic parameters that could be real- 
ized in straightforward adaptions of existing BEC ap- 
paratuses 0, 0|. It was shown that all phenomena 
leading to instability, namely, two-dimensional primary 
collapse, three-dimensional primary collapse, explosion 
of individual solitonic pulses brought about by the lon- 
gitudinally expulsive harmonic trapping potential, and 
secondary collapse caused by soliton-soliton interaction, 
could be avoided by the proper choice of parameters. 
Typical parameters were 10 4 particles, a final scatter- 
ing length of a ~ — 3ao, and trapping frequencies on the 
order of 2tt x 2.2 kHz by 2tt x 2.2 kHz by 2ni x 2.5 Hz. 
After formation via modulational instability seeded by a 
combination of self-interference of the condensate order 
parameter and noise due to the presence of Boguliubov 
quasiparticles and fluctuations in the trapping potential, 
propagating solitonic pulses self-cool to T = on a time 
scale of l/|w z | through the emission of a fraction of a 
percent of the total number of particles 0, 0| . 

In most previous work on attractive Bose-Einstein con- 
densates, reg imes or cycles of runaway instability were 
explored |53t IH^j . Even in the cases where a stable fi- 
nal state was produced, as for example in Refs. [Tsl fli^. 
the majority of the atoms were lost to collapse. In con- 



10 



trast, we have here suggested a way to avoid collapse 
entirely and take advantage of the instabilities inherent 
in switching the interactions in a BEC from repulsive to 
attractive to produce a useful device: namely, a pulsed 
atomic soliton laser. 
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